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ABSTRACT 

Presence  of  surface  imperfections  may  have  significant  impact  on  the  stability  of  the  adjacent  shear  layer  and 
on  the  laminar-turbulent  transition  process  in  this  layer.  The  main  objective  of  the  present  study  is  the 
quantitative  assessment  of  such  effects  in  the  case  of  distributed  surface  roughness.  The  analysis  is  focused  on 
the  linear  stability  theory’  and  considers  asymptotic  as  well  as  transient  growth  of  flow  perturbations.  The 
explicit  calculations  are  carried  out  in  the  context  of  channel  flow. 


1.0  INTRODUCTION 

Understanding  of  the  transition  process  in  flows  over  rough  surfaces  is  one  of  fundamental  questions  in  fluid 
dynamics.  This  question  is  of  practical  interest  in  aerospace  applications  including  the  classical  problem  of  the 
design  of  laminar  airfoils  as  well  as  in  the  more  current  problem  of  the  design  of  surface  modifications  for 
flow  control.  The  latter  application  area  involves  development  of  strategies  for  improvement  of  mixing  for 
better  control/manipulation  of  the  separations  process  in  boundary  layers  and  has  significant  potential  for 
creation  of  means  for  active/passive  manipulation  of  heat  transfer  in  engines. 

Flows  over  rough  walls  have  been  studied  since  the  early  works  of  Hagen  [Hagen  1854]  and  Darcy  [Darcy 
1857],  which  were  focused  on  turbulent  regimes.  Reynolds  [Reynolds  1883]  was  the  first  to  pose  the  problem 
in  the  context  of  laminar-turbulent  transition.  In  spite  of  many  attempts  [Schlichting  1979]  the  resolution  of 
this  problem  is  still  uncertain.  The  original  investigations  involved  measurements  of  flows  in  open  channels 
and  in  pipes.  Various  possible  roughness  forms  were  classified  using  the  concept  of  “equivalent  roughness”. 
Jimenez  [Jimenez  2004]  gives  a  nice  description  of  the  history  and  the  current  state-of-the-art  of  this  concept. 
The  phenomenological  effects  of  the  “equivalent  roughness”  are  summarized  in  the  form  of  friction 
coefficient  [Nikuradse  1933],  [Colebrook  1939],  [Moody  1944].  These  and  other  similar  investigations  show 
that  surface  roughness  contributes  directly  to  the  dynamics  of  flow  only  if  the  wall  is  hydraulically  rough.  The 
concept  of  hydraulic  smoothness  is  conceptually  very  appealing,  however  no  precise  criterion  exists  that 
would  allow  to  predict  if  a  given  surface  can  be  considered  as  being  hydraulically  smooth  for  flow  conditions 
of  interest.  While  the  modelling  concepts  of  this  type  have  been  continuously  re-evaluated  [Bradshaw  2000] 
and  improved  [Waigh  1998],  they  failed  so  far  to  uncover  the  mechanisms  that  govern  the  complex,  flow- 
condition-dependent  interaction  between  the  roughness  geometry  and  the  moving  fluid. 

There  were  a  number  of  concepts  pursued  in  the  attempt  to  link  distributed  surface  roughness  with  the 
transition  process.  Singh  [Sigh  1972],  Lessen  [Lessen  1976]  and  Kendall  [Kendall  1981]  considered 
roughness  -distorted  average  velocity  profile,  but  linear  stability  analysis  of  this  profile  proved  inconclusive. 
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ORGANIZATION 


Merkle  [Merkle  1977]  considered  concepts  based  on  the  roughness-induced  additional  mixing.  Reshotko 
[Reshotko  1984]  and  Morkovin  [Morkovin  1990]  attempted  to  explain  roughness  effects  in  terms  of  Tollmien- 
Schlichting. 

Floryan  (1997)  proposed  a  completely  different  approach  where  roughness  geometry  is  modelled  using 
spectral  expansions;  this  formulation  reduces  the  analysis  of  the  effects  of  different  geometries  to  scans  of 
parameter  space  formed  by  the  coefficients  of  the  expansions.  The  problem  is  split  into  two  steps.  In  the  first 
one  the  basic  state  that  accounts  for  the  spatial  structure  of  the  roughness  is  determined  numerically.  In  the 
second  step,  linear  stability  analysis  that  must  account  for  the  spatial  structure  of  the  roughness  beyond  a 
simple  change  of  the  average  velocity  profile  is  carried  out.  The  above  process  is  illustrated  using  the  simplest 
possible  geometry,  i.e.,  channel  with  one  wavy  wall.  Such  wall  geometry  is  represented  in  terms  of  a  single 
Fourier  mode  and  thus  forms  a  very  convenient  reference  case. 

The  laminar-turbulent  transition  occurs  due  to  nonlinear  effects  induced  by  growing  disturbances.  In  the  case 
of  asymptotic  instability  disturbances  grow  exponentially  without  any  limits  and  always  induce  nonlinear 
effects.  In  the  case  of  transient  growth  disturbances  reach  a  certain  magnitude  but  criterion  that  would  permit 
to  determine  if  this  magnitude  is  sufficient  to  induce  nonlinear  effects  is  not  known.  Both  growth  mechanisms 
will  be  considered  in  this  study. 

This  paper  is  organized  as  follows.  Section  2  provides  brief  description  of  the  basic  state.  Section  3  provides 
formulation  of  the  asymptotic  stability  problem.  Section  4  provides  formulation  of  the  transient  growth 
problem.  Section  5  provides  discussion  of  selected  results.  Section  6  gives  a  short  summary  of  the  main 
conclusions. 


2.0  BASIC  STATE 

The  following  presentation  is  limited  to  a  short  outline.  Details  are  given  by  Floryan  [Floryan  1997].  Consider 
steady,  two-dimensional  flow  in  a  channel  with  corrugated  lower  wall  and  smooth  upper  wall.  The  flow  field 
is  represented  in  the  form 


V2  (x)  =  [u  2  (x,  y),  v2  (x,  y)]  =  V0  (x)  +  Vj  (x) 
P2  (x)  =  Po  (x)  +  Pi  00=  -2x/Re  +  Pi  (x), 


l-y2,0  +  [u1(x,y),v1(x,y)]> 


(2.1) 


where  Vo,  po  describe  the  reference  plane  Poiseuille  flow  and  Vi,pi  describe  modifications  associated  with  the 
presence  of  the  corrugation.  The  channel  extends  from  -go  to  +oo  in  the  x-direction  and  the  locations  of  the 
lower  and  upper  walls,  i.e.,  yL(x)  and  yu(x),  are  defined  as 

yL(x)  =  -l+(Seiax  +  CC),  yu(x)  =  1  ,  (2.2) 

where  CC  stand  for  the  complex  conjugate,  and  S  and  a  denote  the  amplitude  and  the  wave  number  of  the 
wall  corrugation,  respectively.  The  flow  is  scaled  using  the  maximum  of  the  x-velocity  of  the  reference 
Poiseuille  flow  and  the  average  half-channel  height  h.  The  flow  modifications  can  be  represented  using 
Fourier  expansion  in  the  form 

n=+oo 

>f(x,y)=  2>"(y)eln“  (2.3) 

n=— oo 
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where  'P  represents  the  stream  function  defined  in  the  usual  manner,  i.e.,  u ,  =dy'¥,  =  -~SX  T,  On  =  ®*n 

and  star  denotes  complex  conjugate.  Elimination  of  pressure  and  use  of  (2.3)  reduces  the  field  equations  to  the 
equations  for  the  functions  ®n,  n>0,  in  the  form 

D2  -inaRe(u0Dn  -D2u0  ¥  n  -i«Re  X[kDOn-kDk0k  -(n-k)®n-kDkDOk]=0>  (2-5) 

k=-ao 

where  D=d/dy  and  D„=D2-n2a2.  The  boundary  conditions  at  the  channel  walls  are  expressed  in  the  following 
form 

uo(yL(x))+ui(x>yL(x))=°>  vi(x^yL(x))= 0 »  (2-6a) 

u1(x,l)  =  0,  v1(x,l)  =  0.  (2.6b) 

The  required  closing  condition  [Floryan  1997]  has  the  form  of  the  fixed  mass  flux  constraint  and  can  be 
expressed  as 

(y  u  M)  +  ^  y  u  (x))  =  F  +  Q>  (2-7a) 


'F0(-l)  +  'F(x,-l)  =  F,  (2.7b) 

where  T0  =-y3/3+y+2/3  denotes  the  stream  function  of  the  Poiseuille  flow  (continued  analytically  in 
—  l<y<l  +  2S),  Q  stands  for  the  (specified)  volume  flux  and  F  denotes  an  arbitrary  constant  associated  with 
introduction  of  the  stream  function  (F  denotes  value  of  the  stream  function  at  the  lower  wall).  In  the  present 
analysis  F=0  and  Q=4/3.  The  above  problem  is  solved  numerically  with  spectral  accuracy  using  Chebyshev 
expansions.  The  boundary  conditions  at  the  bottom  wall  are  implemented  using  the  immersed  boundary 
conditions  concept  [Szumbarski  1999] [Floryan  2002]  while  the  boundary  conditions  at  the  upper  wall  are 
implemented  using  a  version  of  the  tau  technique. 


3.0  ASYMPTOTIC  LINEAR  STABILITY  PROBLEM 

The  analysis  begins  with  the  governing  equations  in  the  form  of  vorticity  transport.  Details  of  the  formulation 
can  be  found  in  Floryan  [Floryan  1997].  Unsteady,  three-dimensional  disturbances  are  superimposed  on  the 
mean  part  in  the  form 

to  =  co2  (x,  y)  +  co3  (x,  y,  z,  t),  v  =  v2(x,y)  +  v3(x,y,z,t),  (3.1) 

where  subscripts  2  and  3  refer  to  the  mean  flow  and  the  disturbance  field,  respectively.  The  assumed  form  of 
the  flow  is  substituted  into  the  governing  equations,  the  mean  part  is  subtracted  and  the  equations  are 
linearized.  The  resulting  linear  disturbance  equations  have  the  form 

- F  +  (v2- V)co3-(co3- V)v,  +  (v3- V)co,  — (co2- V)v3  =  Re  co3  ,  (3.2a) 

dt 

V-v3=0  ,  co3=Vxv3  (3.2b, c) 
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and  are  subject  to  the  homogeneous  boundary  conditions 
v3(x,l,z)  =  0,  v3(x,yL(x),z)  =  0 

where  yL  is  given  by  Eq.(2.2).  The  mean  flow  has  the  form  (Section  2) 


n=  00.- 

!  (x,  y)  =  [u0  (y),0,o]  +  Yj  [fun)  (y)>  fvn)  (y)>  oJexp(inax) 


(3.2d) 


(3.3) 


where  represent  solution  to  the  problem  (2. 5)-(2. 7)  and  f„n)  =  (fL,  n))*,  f\,n)  =  (f '  n))*. 

Partial  differential  equations  (3.2)  have  coefficients  that  are  functions  of  y  and  x,  being  periodic  in  x.  The  t 
and  z  dependence  can  be  separated  in  the  usual  manner  leading  to  the  solution  in  the  form 


y  3  (x,  y,  z,  t)  = 


rn=+op 

ikm,( 


<y),g<v",(y).gr(y> 


(m)/ 


gi[(8+ma)x+Pz-ot] 


(3.4) 


where  8  and  [1  are  real  and  account  for  the  streamwise  and  spanwise  periodicity  of  the  disturbance  field, 
respectively.  The  exponent  a  is  complex  and  its  imaginary  and  real  parts  describe  the  rate  of  growth  and  the 
frequency  of  the  disturbances,  respectively. 


Substitution  of  (3.3)  and  (3.4)  into  the  disturbance  equations  (3.2)  and  separation  of  Fourier  components 
results  in  a  system  of  ordinary  differential  equations  governing  g<Llm|,g(vni'l,glj’1) .  It  is  convenient  for  analysis 
purposes  to  express  these  relations  in  terms  of  the  y-components  of  disturbance  vorticity 
0<m)  =  tmg^m)  -|3gum)  and  disturbance  velocity  gym),  i.e., 


S(m)0(m)  +  Cg(m)  =  Rg  £  (E  (m,n)g(m-n)  +  E(m,n)0(m-n) }  #  (3.5a) 


T(m)g(m)=_Re^(H 


(m,n) 

V  c 


(m-n) 

V 


H  (m.n)  q  (m— n)  ^ 


-  co  <  m  <  go 


(3.5b) 


where  the  explicit  forms  of  the  operators  E,  El  are  given  in  the  Appendix  A.  The  boundaiy  conditions  have  the 
form 

gf’O)  =  g«(D  =  gf”1  (i)  =  gf"(yL)  =  g<m>(yL)  =  gf"1  (yL)  =  o,  (3,6) 

Conditions  at  the  upper  wall  are  imposed  using  a  variation  of  the  tau  technique,  while  those  at  the  bottom  wall 
are  implemented  using  immersed  boundary  conditions  concept  [Floryan  2002].  Equations  (3.5)  with  boundary 
conditions  (3.6)  have  nontrivial  solutions  only  for  certain  combinations  of  parameters  8,  a  and  p.  The  required 
dispersion  relation  has  to  be  determined  numerically.  For  the  purposes  of  calculations,  the  problem  is  posed  as 
an  eigenvalue  problem  for  a. 
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4.0  LINEAR  TRANSIENT  GROWTH  PROBLEM 

The  analysis  begins  with  the  same  formulation  as  used  in  the  previous  section.  The  following  presentation  is 
limited  to  a  short  outline.  Details  are  given  by  Szumbarski  [Szumbarski  2004].  The  disturbances  are  assumed 
in  the  form 


v3(t,x,y,z)  =  [u(x,y,t),v(x,y,t),w(x,y,t)]el5x+l|3z  =  £  [gSn)(t,y),g$n)(t,y),g^l)(t,y) 


,i(YmX+PZ) 


+  CC  (4.1) 


where  P  and  8  are  real  and  denote  spanwise  and  streamwise  wave  numbers,  respectively,  ym  =  8  +  ma, 
u,v,w  describe  modulation  of  disturbances  by  the  corrugation  and  are  periodic  in  x.  Equations  (4.1)  and  (3.3) 
are  substituted  into  (3.2),  the  unknowns  are  expressed  in  terms  of  the  normal  velocity  component  v3  and  the 
normal  vorticity  component  eoDy,  where 


®Dy  _ 


5U3 

dz 


5w3 

dx 


m=+  oo 

=  -i  X  9(m)(t,y)ei(YmX+Pz), 

m=  -oo 


e(m>=-pgim)+Tmg 


(m) 

w 


(4.2) 


and  the  Fourier  modes  are  separated  resulting  in  the  following  form  of  disturbances  equations  for  0(m),gym), 
me  (-go,  +go) 


(d2 -k2m)dtg^ 


00  / 

Re“'Z(^(vm’n)g(vm+n)  +G(vm’n)gim_n)  +G[)m>n)0(m+n)  +G(nn,n,0(m-n) 
n=l 


)+G<vo)gr 


(4.3a) 


5t0(m)  -  Re-1  Id2  -k^  -iReymu0  0(m)  +  pDu0  g<m)  = 


(4.3b) 


=  Re 


1 £  (sf ’n)gf+n)  +  sim-n)gtm“n)  +  Sgm,n,0(m+n)  +  s<m’n)0(m-n) )+  sfgtm)  +  SgO)0(m)  . 


n=l 


The  explicit  forms  of  the  operators  Gv,G0,Gv,G9,Sv,Sg,Sv,S9are  given  in  the  Appendix  B.  It  can  be 
shown  that  boundary  conditions  (3.3)  take  the  form 


gum)  (t,  yL  00)  =  o,  gtm>  (t,  yL  (x))  =  0,  g^  (t,  yL  (x))  =  0, 0(m)  (t,  yL  (x))  =  0 

g^m)(t,i)  =  o,  glm)(t,i)  =  o,  g^CU)  =  0, 0 (m) (t, l )  =  o . 


(4.4a) 

(4.4b) 


where  me(-Go,+Go).  The  boundary  conditions  at  the  upper  wall  are  implemented  using  a  version  of  the  tau 
technique,  while  boundary  conditions  at  the  lower  wall  are  implemented  using  the  immersed  boundary 
conditions  concept.  Problem  (4.3)-(4.4)  has  to  be  supplemented  by  the  initial  conditions.  The  solution 
involves  search  for  the  form  of  the  initial  conditions  that  gives  the  largest  possible  transient  growth  (optimal 
disturbances). 
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5.0  DISCUSSION  OF  RESULTS 

The  unstable  disturbances  have  the  form  of  either  streamwise  vortices  or  travelling  waves  [Floryan  2003].  The 
vortices  are  induced  by  the  presence  of  the  roughness,  but  only  when  the  roughness  amplitude  and  wave 
number  are  in  the  correct  range.  Figure  1  displays  a  typical  stability  diagram  for  the  streamwise  vortices. 
Reduction  of  either  the  roughness  amplitude  S  or  the  flow  Reynolds  number  Re  reduces  the  unstable  zone. 
Critical  conditions  are  found  when  this  zone  shrinks  to  a  single  point. 


Figure  1.  Amplification  rates  lmag(o)  of  disturbances  in  the  form  of  streamwise  vortices  as  a  function  of  the 
roughness  wave  number  a  and  the  vortex  wave  number  p  for  the  flow  Reynolds  number  Re=5000  and  the 
roughness  amplitude  S=0.0075. 

Travelling  waves  become  unstable  when  the  Reynolds  number  reaches  the  critical  value  of  Re=5772  in  a 
smooth  channel.  Presence  of  roughness  significantly  reduces  this  value.  Figure  2  displays  a  typical  stability 
diagram  in  the  case  of  rough  channel.  The  size  of  the  unstable  zone  decreases  when  either  the  roughness 
amplitude  S  or  the  flow  Reynolds  number  Re  decrease.  The  critical  conditions  are  found  when  this  zobe 
shrinks  to  zero. 


Figure  2.  Curves  of  constant  amplification  lmag(o)  for  disturbances  in  the  form  of  two-dimensional 

travelling  waves  for  the  roughness  amplitude  S=0.0085  and  the  flow  Reynolds  number  Re=5600  as  a  function  of 
the  roughness  wave  number  a  and  the  disturbance  wave  number  S.  Dash  lines  correspond  to  Real(o). 

Figure  3  displays  critical  curves  for  both  types  of  instability.  The  area  below  these  curves  defines  conditions 
when  the  wall  behaves  as  hydraulically  smooth.  There  will  be  at  least  one  roughness  wave  number  that  would 
induce  instability  for  a  given  roughness  amplitude  when  Reynolds  number  increases  above  the  critical  value 
defined  by  these  curves. 
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Figure  3.  Variations  of  the  global  critical  Reynolds  number  Reg,cr  describing  the  travelling  wave  instability  and 
the  vortex-like  instability  as  a  function  of  the  roughness  amplitude  S.  The  area  below  and  to  the  left  of  both  curves 
corresponds  to  the  flow  conditions  and  the  wall  geometry  that  do  not  produce  any  instability.  Such  wall  operating 
under  the  specified  Reynolds  number  behaves  as  the  hydraulically  smooth  wall. 

The  definition  of  the  hydraulically  smooth  wall  stated  above  is  based  on  the  asymptotic  instability.  Such 
instability  always  induces  transition.  The  reader  should  note  that  this  definition  omits  transient  growth  that 
may  induce  transition  under  subcritical  conditions.  The  formulation  presented  in  Section  4  provides  means  for 
assessing  the  strength  of  the  transient  growth.  This  issue  is  currently  subject  to  detailed  investigation. 

6.0  SUMMARY 

Formulation  of  the  problem  of  growth  of  linear  disturbances  in  a  rough  channel  is  presented.  The  roughness 
geometry  is  represented  in  terms  of  Fourier  series.  Two  types  of  disturbance  growth  are  considered,  i.e., 
exponential  growth  associated  with  the  asymptotic  instability  and  transient  growth.  Detailed  results  are  given 
in  the  case  of  roughness  described  by  a  single  Fourier  mode  and  asymptotic  instability.  Conditions  under 
which  roughness  is  not  hydraulically  active  have  been  identified  using  criterion  based  on  the  roughness 
ability  to  induce  asymptotic  instability. 

7.0  APPENDIX  A 

Operators  used  in  Eqs.  (3.6)  have  the  following  definitions: 

E(m,n)  =p(-Df(n)  +inak-2_nf(n)D2);  =  itmf <n>  +  (1  +  natm_nk-2_n)ff  D  ; 

=inak-ln(p2  -tmtm_n)Dff  D  +  k^k-2_n(p2  +  tm_ntm„2n )ff  D  + 
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ORGANIZATION 


8.0  APPENDIX  B 

Operators  used  in  Eqs  (4.3)  are  defined  as  follows: 
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